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Abstract 

The spectrum of the Dirac Hamiltonian in the background of crossing vortices 
is studied. To exploit the index theorem, and in analogy to the lattice the space- 
time manifold is chosen to be the four-torus T 4 . For sake of simplicity we consider 
two idealized cases: infinitely fat and thin transversally intersecting vortices. The 
time-dependent spectrum of the Dirac Hamiltonian is calculated and in particular 
the influence of the vortex crossing on the quark spectrum is investigated. For the 
infinitely fat intersecting vortices it is found that zero modes of the four-dimensional 
Dirac operator can be expressed in terms of those eigenspinors of the Euclidean time- 
dependent Dirac Hamiltonian, which cross zero energy. For thin intersecting vortices 
the time gradient of the spectral flow of the Dirac Hamiltonian is steepest at the time 
at which the vortices cross each other. 
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1 Introduction 



Center vortices offer an appealing picture of confinement PjJ |2] , for a recent review see [3] . 
Furthermore, the deconfinement phase transition can be understood as a depercolation 
transition jlj and the free energy of a center vortex vanishes in the confinement regime 
[Hj. In fact there is mounting evidence that center vortices play an important role in the 
QCD vacuum. Lattice calculations show that when center vortices are removed from the 
Yang-Mills ensemble the string tension |6 b as well as the quark condensate and topolog- 
ical charge [7j are lost 1 . This indicates that confinement, spontaneous breaking of chiral 
symmetry and the chiral anomaly are all triggered by center vortex configurations in the 
Yang-Mills vacuum. Furthermore, lattice calculations [B] show that restoration of chiral 
symmetry occurs at finite temperature exactly at the deconfinement phase transition. We 
can therefore expect that confinement and spontaneous breaking of chiral symmetry are 
related phenomena. While the emergence of the confining potential, i.e. the string ten- 
sion, as well as the appearance of the deconfinement phase transition can be rather easily 
understood in a random center vortex model in an analytic fashion, the mechanism of spon- 
taneous breaking of chiral symmetry has not been explained in the vortex picture of the 
QCD vacuum, yet. The general folklore is that spontaneous breaking of chiral symmetry is 
generated by the quasi-zero modes in topologically non-trivial field configurations, which 
give rise to a quasi continuous Dirac spectrum with non-zero level density p(0) at zero vir- 
tually and, by the Banks Casher relation, to a non-zero quark condensate < qq >~ p(0). 
In fact these calculations [9 j show, that a non-zero level density at zero virtuality occurs in 
the confinement phase while a finite gap in the Dirac spectrum is found above the critical 
temperature. 

It is therefore clear that the topological properties of center vortices and their manifes- 
tation in the quark spectrum will play the key role in the understanding of spontaneous 
breaking of chiral symmetry in the center vortex picture of the QCD vacuum. In D = 4 
space-time dimensions center vortices represent closed surfaces of quantized electromag- 
netic flux which produce a non-trivial center element for a Wilson loop when they are 
topologically non-trivially linked to the latter The topological charge of center vor- 
tices is given by their (oriented) intersection number IT3*1 ITlj. According to the index 
theorem a non-zero topological charge implies zero modes of the Dirac operator |T5]. In 
[To] the quark zero modes in transversally intersecting vortex surfaces have been studied. 
It was found that these modes are localized at the vortex sheets, in particular at their 
intersection points. In D = 3 dimensions a center vortex represents generically a time- 
dependent closed flux loop. In [T2] it was shown that the topological charge of generic 
center vortices can be expressed by the temporal change of their writhing number, and 
furthermore transversal intersection points in D = 4 correspond to crossing of vortex lines 
in D = 3. Thus vortex crossings carry topological charge and by the index theorem should 
consequently manifest themselves in the quark spectrum. 

In this paper we study the quark spectrum in the background of time-dependent crossing 

1 In fact the finding of the lattice calculation concerning string tension and deconfinement phase tran- 
sition can be reproduced in a center vortex model [TO1 ITTj . 
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center vortices. In particular we investigate how the crossing of vortices manifests itself 
in the quark spectrum. This might bring some light on the mechanism of spontaneous 
breaking of chiral symmetry in the vortex picture, since it is believed that chiral symmetry 
breaking is tightly related to the topological properties of gauge fields. 

Throughout this paper the space-time- manifold is chosen to be the four-torus T 4 . There 
is a variety of reasons for studying T 4 : first, T 4 allows to use the Atiyah-Singer index 
theorem in a stringent fashion and allows for normalizable spinors in the background of 
integer flux. Second, the torus simulates a periodic arrangement of vortices. This is much 
closer to a percolated vortex cluster than two isolated vortices in M 4 . Finally, the torus is 
the space-time manifold that is also used in lattice calculations. 

The quark spectrum and its spectral flow in time-dependent background gauge fields 
have been addressed in several publications before, e.g. [13 CHI HE]- In these papers strict 
relations between the spectral flow of the Dirac Hamiltonian and topological properties 
of the corresponding background gauge field have been derived in Minkowski space-time. 
In [T71 UH] the spectral flow was considered on an infinite time interval, where the back- 
ground gauge field has to become a pure gauge at t — — oo and t = +00. In [19 j only 
spherically symmetric gauge fields have been examined, but the derived results are also 
valid on a finite time interval. In this paper we analyze the spectral flow in the presence 
of intersecting center vortices on a compact Euclidean space-time manifold. The present 
paper is a continuation of ref. JH] where the quark zero modes in intersecting center vortex 
sheets were considered. In the present paper we will not confine ourselves to the quark 
zero modes but study the spectral flow as function of the Euclidean time for the time- 
dependent vortex crossings. From studying the spectral flow we hope to get some insight 
into the mechanism which produces the quasi-continuous spectrum of Dirac states near 
zero virtuality. As a byproduct we will also establish a strict relation between the quark 
zero modes and the (time-dependent) eigenfunctions of the Dirac Hamiltonian in the case 
where the background fields consist of "vortices" with homogeneously spread out magnetic 
flux. 

The organization of the paper is as follows: In the next section we consider the math- 
ematical idealization of infinitely thick transversally intersecting vortices, i.e. we consider 
"vortices" which are homogeneously distributed over the transversal dimensions, so that 
they represent homogeneous field configurations of constant electromagnetic flux. We con- 
sider two orthogonally homogeneously distributed fluxes, which can be considered as two 
infinitely thick intersecting vortices, where the intersection "point" is smeared out over 
the whole universe. Accordingly the topological charge is homogeneously distributed over 
space-time. We will study then the eigenmodes of the Dirac Hamiltonian as function of 
the adiabatic time-evolution of the vortex flux assuming that the three-space is given by 
a three-torus. In section |3] we study the opposite case of thin transversally intersecting 
vortices assuming that the two intersecting vortices are orthogonal to each other. We study 
the spectral flow of the quarks in this vortex background. Finally in section 0] we allow 
the intersecting vortex sheets to be non-orthogonal and calculate again the spectral flow 
of the quarks. Some concluding remarks are given in section |S| 
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2 Fermions in homogeneously intersecting "vortex" 
fields 

The topological charge of a gauge field configuration 

v = / trFF (2.1) 

32tt 2 J v ; 

consisting of center vortices is located at the intersection points of the vortices [T3"| IT2*] . In 
[T2*] the evolution in time of the topological charge of crossing center vortex configuration 
has been studied. It is the aim of this paper to investigate the spectral flow of the Dirac 
Hamiltonian in the background of crossing vortices. 

As space-time manifold we choose the four-torus T 4 . It can be seen as M 4 modulo the 
lattice 

x G 1R 4 | x = n^e^ , G Z 

where {e M } , /i = 1, . . . ,4 are four orthogonal basis vectors with lengths — |c/x|- 

To demonstrate how the spectral flow emerges for intersecting center vortices we start 
with the following SU (2) gauge field configuration on the four-torus T 4 

A 1 = -2mx 2 /(L x L 2 )T Z , A 2 = , A 3 = -2mx 4 /{L 3 L 4 )r 3 , A A = . (2.2) 

This gauge potential describes two intersecting "vortices" whose fluxes are spread out over 
the whole x x -x 2 - and x 3 -x 4 -, respectively, planes. Indeed, this configuration has constant 
non-zero field strength components 

= 27ri/(L 1 L 2 )T 3 , F 34 = 2m/{L 3 L A )r 3 (2.3) 

so that their fluxes are given by 



/ 
/ 



dS 12 F 12 = I 1 dxiA!(x 2 = L 2 ) = -%2itt 3 , (2.4) 
Jo 

dE 34 F34 = / dx 3 A 3 (x 4 = L 4 ) = -i2v:t 3 , (2.5) 







(2.6) 



which is twice the flux of a center vortex. The line integrals appearing here represent the 
only non- vanishing contributions from the Wilson loops encircling the whole X\-Xt and x 3 - 
x 4 -plane, respectively, see fig. [TJ We choose the flux of these configurations to be twice the 
flux of a center vortex since a single center vortex cannot be accommodated by the torus 
with fermions in the fundamental representation of the gauge group, see the discussion 
given below eq. (|2.9|) . 

Furthermore, the configuration (j2.2j) carries topological charge v = —2 which is homo- 
geneously distributed over space-time, i.e. the intersection "point" of these two very fat 
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u x 1 

Figure 1: Wilson loop encircling the 
whole Xi-x 2 -plane. 




Figure 2: Illustration of the spectral flow 
of left and right handed fermions in topolog- 
ically non-trivial constant background field 
strength as function of time r. 



vortices is smeared out over the whole space-time. The gauge potential ()2.2|) satisfies the 
boundary conditions 

A^x + eJ = A u u »{x) = U^x)A v {x)U;\x) + U^x)d v U;\x), (2.7) 

where the transition functions can be chosen as follows: 

Ui = U 3 = t , U 2 = exp (2irixi/L 1 T 3 ) , U A = exp (2111x3/ L 3 t 3 ) ■ (2.8) 

In general the transition functions have to fulfill the cocycle condition 

U v (x + e^)U^x)Z^ v = U^x + e v )U v {x) , (2.9) 

where are elements from the center of the gauge group. Throughout this paper we 
consider fermions in the fundamental representation of SU(2). This implies that all 
are equal to 1, because fermions in the fundamental representation cannot tolerate twisted 
boundary conditions which imply the presence of a single center vortex. For later use we 
also note that under a gauge transformation — > A^ the transition functions transform 

U%(x) = V(x + e^U^Vix)- 1 . (2.10) 

Furthermore, since our gauge configuration ()2.2j) satisfies the Weyl gauge A 4 = its topo- 
logical charge (|2.1j) can be expressed in terms of the (in general non-Abelian) Chern-Simons 
action Scs(x4> as [12] 

u= J * dx^Scsix^) . (2.11) 
J 

For well localized center vortices (representing closed flux loops in D = 3) the Chern-Simons 
action is nothing but their writhing number ^2]- For the particular gauge potential (j2.2j) 
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the Chern-Simons action reduces to 

Scsix,) = ~ [ tr(A 3 F 12 )d 3 x = -2^i := -2r (2.12) 

so that eq. (|2.11|) yields indeed v = —2 when the vortex evolves from r = to r = 1. 

We are interested in the time-evolution of fermions in the intersecting center vortex 
background described by the gauge potential (|2.2)l . Specifically we want to study the 
time-evolution of the spectrum of the three dimensional Dirac-Hamiltonian 

H(r)ip = E(T)ip, 

i.e. the Dirac energies E(t) as function of the dimensionless time r. 
The Dirac Hamiltonian of a massless fermion reads 2 

tf = - 7l7iDi =(-7 D < J D _), (2.13) 

where Dj = d{ + Ai is the covariant derivative. Under gauge transformations A — > A v the 
Dirac spinors in the fundamental representation of 577(2) transform as 

ij){x) -> ip v (x) = V(x)ip(x) . (2.14) 

In view of equation (J2.7|) this implies the following boundary conditions for fermionic fields 

+ e M ) = U^x)^{x) . (2.15) 

Because the gauge potential (J2.2)) and the transition functions (J2.8)) are purely Abelian 
the two color components of the fermionic field decouple from each other. Therefore it is 
sufficient to consider only one of the two color components, i.e. we restrict ourselves to 
spinors with the lower color component equal zero. To get the desired results for spinors 
with the opposite color one has to change the sign of the energy and to exchange the upper 
with the lower color component. If not stated otherwise, in the following we will always 
consider the upper color component of the Dirac field. 

Coming back to the Dirac spectrum one first observes that it is invariant under the 
shift r — > r + 1. This is because the two corresponding gauge potentials are related by the 
periodic gauge transformation V = exp (— 2nix 3 / L 3 t 3 ) which leaves the above introduced 
transition functions Uy^/s invariant, see eq. (|2.1(J|) . 

Since H commutes with 75 the eigenvectors of H can be chosen as eigenvectors of 75. 
Hence the eigenvectors of H carry good chirality. It remains to find the eigenvectors ip of 
the Dirac operator 

p (3) =i(T i D i ,p i3) $ = E$ (2.16) 

in three dimensions. Because of the periodicity in the 3-direction (U 3 = 1) one can make 
the ansatz 

r tp n (xi,x 2 ,x 3 ) = ip n (xi, x 2 ) exp (-2irinx 3 / L 3 ) . 

2 The representation of the Dirac matrices used in this paper can be found in Appendix 1X1 
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Inserting this ansatz into the eigenvalue equation (|2.1fij) the latter reduces to the eigenvalue 
problem 

[p (2) + E n a 3 j ij n = A^n , E n = 2vr-^— , (2.17) 

where p^ = i (o^Dx + o" 2 D 2 ) is the massless Dirac operator in two dimensions. For the 
gauge potential (|2.2|) the eigenvectors and eigenvalues of p^ are known [20 : 



P {2) <f) m = E m (f) m , E m = sign(m)W j-j- , m E Z. (2.18) 

According to the index theorem for Abelian gauge fields in D = 2 dimensions, (for each 
color component) there is exactly one zero mode of P^, namely 4> m =o- This zero mode is 
left handed, i.e. only the lower component of this two-spinor is non-zero. Therefore, the 
two-spinor ^) nm= Q = (p m= o(xi, x 2 ) exp (— 2ninx 3 /L 3 ) is eigenvector of p^ with eigenvalue 

— E n , see eq. f!2.17|) . The explicit form of the function 4> m= o is written down in eq. (jB.4|) . 
For m/fl the eigenvectors ip n ()2.17|) are given by linear combinations of 4> m and 0_ m 

Putting this ansatz into eq. ()2.17|) and noting that 3 cr 3 (j) m = 0_ m one gets the following 
eigenvalue problem for the vector (a nm , f3 nm ): 

Em, E n \ I Ot nm \ I Ot nm \ /o 1 n\ 

Enm a • (2.19; 



E n E m J V /3 nm J \ P; 

This eigenvalue problem has for fixed (n, m) two solutions: 



E nm = ±yjE 2 m + E 2 , (2.20) 
where the sign is given by the sign of E m . One of the two solutions corresponds to positive 



m E 7L and the other to negative m E Z. Therefore we set E nm = E m \J 1 + (E n /E m ) 2 , m ^ 
0. The corresponding coefficients a nm and (3 nm can be calculated from ()2.19|) . 

Using the tensor product notation 4> ® Xi f° r right handed {I = 1) and left handed 
(1 = 2) four-spinors 

= ^ j , 0®X2= y j , (2.21) 
the eigenvalues and eigenfunctions of the Dirac Hamiltonian H ()2.13j) can be written as 



4n 2 (n + t) 2 

E mnl = { V L l E2 m ) , (2.22) 

— \n + t) , m = 

L 3 



This is because p( 2 ) and 03 anti-commute. 
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Figure 3: Probability density of the mode ^002 at con- 
stant x 3 . 



^Pnrrd 



0o ® Xi e { - 2mnx ' i,Li) 




j <8> Xi e { ~ 2mnX3/Ls) 
, n e Z, I = 1,2. 



, m,neZ,m/0,i = l 1 2 l (2.23) 

(2.24) 



Obviously, only the m = level can cross zero energy. As the time r evolves from r = to 
t = 1 the left handed (/ = 2, m = 0) modes show a spectral flow from positive to negative 
energy and the right handed (I = 1, m = 0) modes from negative to positive energy, and 
there is a level crossing through zero energy for n = — r, see fig. 121 The eigenmodes of the 
lower color component show the opposite behavior. 

Semi-classically, one can interpret the spectral flow induced by the time-evolution r — > 
r + 1 as the creation of a (right handed) particle and a (left handed) anti-particle, i.e. the 
number of right-handed particles increases by 2. Correspondingly the Chern-Simons action 
S cs (r) changes by -2. 

The probability density of the spinor ip n=0tm=0 j =2 is plotted in fig. El for fixed x 3 . This 
density is independent of £3 and of the adiabatic parameter r. The probability density has 
a zero and, obviously, breaks the translational invariance of the field strength 4 (J2.3j) . In 
P^] it has been shown that the position of this zero can be moved by adding a constant 
background gauge potential which does not change the field strength and the boundary 
conditions. 

Let us now come back to the Dirac operator on the four-torus. From the index theorem 
in D = 4 it is known that there is exactly one left handed zero mode ipo of the Dirac 
operator with the above defined gauge potential. This zero mode can be explicitly written 
down JH] and is given in appendix El One can show that this zero mode ipo can be written 
as the sum over all left handed eigenmodes ipnmi of the Dirac Hamiltonian H which cross 
zero energy as the time r passes from r = —00 to r = +00, which are the states with the 

4 This phenomenon is already known from the localized Landau orbits in a constant magnetic field. 
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quantum numbers m — 0, 1 — 2. The time dependence of a solution £ n0 « = fnoi{x4)ip n oi of 
the (Euclidean) time-dependent Dirac equation 



is given by 



d^nOl — H£,n0l (2.25) 



d4,fn0l(x4> = EnOifnOi^Xi) , (2.26) 

r / \ 2 ' 



fnoi(x 4 ) = C^exp^-(-l)V-^(n+-^) ). (2.27) 



The right-handed (/ = 1) solutions are exponentially growing for x 4 — > ±00 whereas the 
left-handed solutions (/ = 2) are exponentially decreasing. Therefore one can construct 
a solution which is quasi-periodic 5 in X4 from the left-handed solutions £„02 by a simple 
summation over all modes n 6 Z, putting C n o2 = l,Vn G Z: 

"T^ -^3 V -^4/ -^3 

^ 3 4 ' n£Z ^ 3 ^ 3 3 ' ' 

= to.^®X 2 «p{-^-je(^ + i^,<^). (2.28) 

34_/ \ 3 3 3/ 

Here (2, z£) is the theta function which is defined in eq. (jA.5|) Comparing with eqs. flB.4IB.5|) 
one observes that this is precisely the zero mode ipo of the Dirac operator on the four-torus. 
In the limit I/4/L3 ^> 1 and la^/L^J <C 1 the part of the sum in eq. (|2.28|) with n = 
dominates, i.e. the sum can be approximated by 1. Neglecting the £4 dependence one 
obtains the eigenvector ^002 which crosses zero energy for x 4 = 0, i.e. the zero mode on the 
four-torus approximates the time-evolution of that adiabatic eigenvector of H who crosses 
zero energy during the evolution. One reason for getting the strict relation (J2.28)) between 
the zero mode on T 4 and the eigenmodes of the Dirac-Hamiltonian on T 3 is that the eigen- 
states ipn02 are independent of x 4 . Therefore the adiabatic approximation becomes exact 
and the X4 dependent adiabatic solutions of the stationary Dirac equation are exact zero 
modes of the Dirac equation on T 3 x R. 



3 Spectral flow for orthogonally intersecting vortices 

In this chapter we consider two thin but smeared out vortices carrying both twice the flux 
of a center vortex in SU(2). We choose one vortex to be parallel to the X3-X4 plane and 

J Thc solution will be periodic under X4 — > £4 + L4 up to a gauge transformation with the transition 
function C/4, see eq. (12. 8|) . 
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located at x\ = xf\x2 = x% and the other parallel to the X1-X2 plane and located at 
4 . The gauge potential of these two vortices can be written down in 



terms of theta functions jTB] . To this end one introduces complex variables 

u = (xi + 1x2) / L 1 , = (xf' +1X2 

v = (23 + i2 4 ) /L 3 , = (xf^+ix^ 



U, V. 




(3.1) 



In terms of the new variable u the gauge potential corresponding to the first vortex 
reads [TH] 

A? = -id 2 (f>(uM \r 12 ) , 4 1} = id 1( j ) (uM°\r 12 ) , A™ = A™ = 0, (3.2) 
where = L u /L^ and the profile function <fr is defined in terms of theta functions 6 : 

cf>(u, « (0) ,r 12 ) = ^ log (e+(u, r 12 )9+(u, uM,m) + 0-(u, u^°\ m)9-{u, r 12 ) 
± (u,u< ( Vi 2 ) = 0(M + ^ + V 12 - M (o) ±e,ir 12 ), £6l + . (3.3) 

The parameter e is a measure for the thickness of the vortex, e 2 is approximately the 
ratio between the transversal vortex extension and the area of the two torus. The gauge 
potential A^ of the second vortex is analogously given by 

4 2) = -id^{v, «(°), r 34 ) , = ids<l>(v, v<®, r 34 ) , A? = A { 2 2) = . (3.4) 

Obviously the gauge potential of the complete vortex configuration is given by the sum 
A = A^ + A^ . For the subsequent considerations it is important that the gauge potential 
of the first vortex is given by components A\ and A2 which are independent of 23, 24 and 
the gauge potential of the second one by components A% and A± which are independent of 
Xi,X2- This means that the covariant derivatives D x and D 2 depend only on 21, 2 2 and the 
covariant derivative D 3 depends only on 23 and 24. Under this assumption the spectrum of 
the Dirac operator D^n = ia^Dk in D = 3 dimensions can be written down if the spectrum 
of the Dirac operator Pi 2 \ = * + ^2^2) on the two torus (spanned by the basis vectors 

e±, e 2 ) is known. 

One further remark is in order here. In Euclidean space-time the Dirac Hamiltonian 
()2.13|) is usually calculated in the Weyl gauge AY = 0. It is possible to write down a gauge 
transformation V which transforms the gauge potential given in eq. (J3.4|) into the Weyl 
gauge: 

(rt=X4 \ 
/ A* (2, t)dt\ , 2 = (21, 2 2 , 23) , (3.5) 

where V denotes path ordering. The gauge transformation V will not change the transition 
functions C/1/2/3, because A4 and therefore also V are Abelian and periodic in these three 

6 Our conventions concerning theta functions can be found in appendix 1X1 
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directions, cf. eq. (j2.1()jl . The gauge transformation V will in general alter the transition 
function U4. But the transformation of U4 is irrelevant in this context, because we are only 
interested in the spectrum of the Dirac Hamiltonian at a constant time x±. Therefore it 
makes no difference to calculate the eigenspinors and eigenvalues of H in the original gauge 
given in eqs. (j3.2H3.4j) or in the Weyl gauge - the difference is simply a gauge transformation 
which is periodic in the spatial directions. 

As in the case of constant field strength one can construct the eigenstates of H with 
the help of the eigenstates {(fi m , m G Z} of the Dirac operator p, 2 % = i (eriDi + O2D2) with 
vortex background gauge potential on the two torus 

P{2)4 > m = E m (f) m . 

As before, from the index theorem in D = 2 dimensions it follows that there is exactly one 
left handed zero mode for the upper color index of D^ on the two torus, i.e. E m = 0, if 

and only if m = (cf. eq. (j2.18jl ). Assuming that the eigenstates {<p m } of D^ are known 
one can make the ansatz (cf. eq. (j2.23l2.24j0 

?Pnml = ((Xnm&m + Pnm4>-m) ® Xl K , m,U G Z, 171 ^ , / = 1,2 , (3.6) 

tpnoi = a n0 (}) ® Xi h n , n G Z, I = 1,2, (3.7) 
where h n is an eigenfunction of 

iD 3 h n = E n h n 

which depends on X3 and £4. Noting the periodicity of h n (because U3 = 1) one obtains 
the solution 

h n (x 3 ,x 4 ) = exp l-2-Kin— + — J dx' 3 A 3 (x' 3 , x 4 ) - J dx' 3 A 3 (x' 3 , x A )\ , (3.8) 
2irn i f L3 

E n = — + — dx' 3 A 3 (x' 3 ) . (3.9) 
Inserting this ansatz into the stationary Dirac equation 

Hlpnml E nm iljj nrn i 

yields (as in the case of constant field strength considered before) the eigenvalue equation 
(I2.19j) for the eigenvectors (a nm , (3 nm ) and the energies E nm [. The energy eigenvalues are 
given by 



E nml = E m ^l + {b/E m )\ m/0 I (310) 

y —b , m = J 

b = 27T— + — / dx 3 A 3 . (3.11) 

-^3 -^3 Jo 
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v (0) Y (l) 
X 4 X 4 



Figure 4: Integration path on the two 
torus. 







x : time at which 
the vortices cross 
\ each other 




X = X4 / L4 




1 1 3* 

x 1 



Figure 5: Spectral flow of the energy 
levels in the case of orthogonally in- 
tersecting vortices. Shown is the left 
(right) handed energy level for the up- 
per (lower) color component. 



Obviously, only eigenvectors with m = can cross zero energy. The integral in eq. ()3.11|) 
is directly related to the field strength component F 34 : 

L s ... r L 3 r L 3 r x ^ 



?>) = / dx' 3 A 3 (x' 3 ,xP 

JO JO Jx 



dx' 3 A 3 (x' 3 ,x\ L ') = I dx' 3 A 3 (x' 3 ,x^) + I dx' 3 I (Q) dx' A F, 3 (x' 3 ,x' 4 ).(3.12) 

4 



Therefore, one obtains 

L3 x m 

E n0 i(x 4 1] ) = E n01 (x 4 0) )-(-l) 1 ^- / 3 d^ H dx' 4 F 43 (x' 3 ,x' 4 ), (3.13) 

L 3 Jo Jx { ° ] 

i.e. the energy E n oi(x 4 ) changes most rapidly where the field strength is peaked and this 
is at the location of the vortex in the x^-x 4 plane (which is of course the time where the 
two orthogonal vortices intersect on T 4 ), see figs. 0] and El It turns out that left handed 
(/ = 2) and right handed (1 = 1) spinors cross zero energy level starting from positive to 
negative energies and vice versa. The time x 4 ^ for which the crossing through zero energy 
takes place can be shifted if one adds a constant to the gauge potential A 3 . This addition 
of a constant will not change the field strength, but it can obviously change the time at 
which the zero energy crossing occurs. Therefore, the precise time x 4 of the level crossing 
seems to be irrelevant. 

Obviously, in the case of crossing vortices the adiabatic approximation is (especially 
in the neighborhood of the vortex) not very good, because the eigenstates ip nm i depend 
strongly on x 4 . Therefore in this case we do not have such a simple relation between the 
zero mode on the four-torus and the adiabatic eigenstates of the Dirac Hamiltonian, see 
eq. (EI2HD - 

The probability density of the eigenvectors ip nm 2 is independent of £3 and of the adiabatic 
parameter x 4 . In fig. El the probability density of the mode 1P002 has been plotted for 
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L2 



Figure 6: Probability density of the 
eigenvector ip 00 2 at constant 2 3 . 



Figure 7: Spectral flow for M = 2 vor- 
tices and N = 3 anti-vortices respec- 
tively located at times = X4//L4. 



constant 23. As in the case of constant field strength there is a zero of the density. But, 
in addition, here the density is peaked at the position of the vortex in the x\ — 22 plane. 

It is easy to generalize the previous results to background fields consisting of more than 
two intersecting vortices. In what follows we consider one vortex parallel to the 23-24 plane 
and located at x\ = xf \ 22 = 22 an d a number M of vortices and N anti- vortices which 
are all parallel to the 2i-2 2 plane, located at 2 3 = 23 , 24 = x£\ k — 1, . . . , (M + N), and 
carrying magnetic flux (p^ = ±2ir. The topological charge of this configuration is given by 
the sum v = j- <pk = M — N over the magnetic fluxes of the individual vortices piercing 
the 2 3 -2 4 plane. The gauge potential corresponding to these vortices is simply given by the 
sum over the gauge potentials of the individual vortices, see eq. (13. 4|) . The eigenf unctions 
and eigenvalues of the Dirac Hamiltonian are given by eqs. ()3. 6113 .11]) with gauge potential 
A3 replaced by the corresponding expression of the multiple vortex configuration. The 
spectral flow is shown in fig. for M = 2 vortices and N = 3 anti- vortices. Obviously, 
the gradient of the spectral flow is maximal at times 24 where two vortices cross each 
other. The sign of the gradient depends on the flux of the particular vortex located at the 
time x£\ cf. eq. ()3.13|) . For a gas of orthogonally intersecting vortices the spectrum of the 
Dirac Hamiltonian fluctuates with time r, cf. eq. (|3.1()jl . For m = (E m = 0) the energy 
eigenvalue fluctuates around zero. 

Let us now see how these fluctuations in the spectral flow manifest themselves in the 
spectrum of the D = 4 Dirac operator. Considering a multi-vortex configuration as given 
above with M = N, i.e. with topological charge v — 0. A numerical investigation shows 
that the module of the smallest eigenvalue of the D = 4 Dirac operator decreases with 
increasing number of vortex-anti- vortex pairs. It seams that the fluctuations of the energy 
eigenvalues of the Dirac Hamiltonian lower the eigenvalue of the Dirac operator with small- 
est absolute value. It remains to be seen whether this observation can provide a mechanism 
to generate a non-zero level density at zero virtuality which is sufficient for spontaneous 
breaking of chiral symmetry. This will be subject to further research. 
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X 3 1-3 A 3 

Figure 8: Location of the two vortices 
on the three-torus. 



4 Spectral flow for non-orthogonally intersecting vor- 
tices 

As before we will consider a configuration on T 4 which describes two intersecting center 
vortices. On a spatial torus T 3 one vortex is static and parallel to the X3-axis located 
at x\ = xf\x2 = x^ . The other vortex is time-dependent and parallel to the rci-axis 
and located at x 3 = xf\x2 = (L 2 /-L 4 )x 4 , see fig. |H1 The two vortices intersect in T 4 at 
= = 1, . . . , 3, x 4 = (L4/ L 2 )xf\ The gauge potential of the second vortex can be 

constructed from the potential of a vortex which is originally parallel to the X\ — X2 plane 
by rotation in the a; 2 -:E 4 -plane. Note that the rotated vortex has magnetic and electric field 
strength (the original vortex only has an electric field strength). To write down the gauge 
potential of the vortex configuration one again introduces complex variables u, v. 

u = (x 1 + ix 2 )/L 1 , m (0) = (4 0) + ix 2 )S ) l L i , / 4 -Q 

v = (x 3 + i (x 4 - x 2 r 24 )) /L 3 , = xf ] / L 3 . 

The gauge potential of the first vortex is again given by ^E] 

A { ? = -td 2 <p{u, u {0 \ r 12 ) , A { 2 1) = td 1( p{u, u (0) , r 12 ) , A { 3 1] = A^ = . (4.2) 

The gauge potential A^ of the second (time-dependent) vortex is slightly more involved: 

A 1 = , A 2 = -iT 2A d 3 <t>{v, t/°>, r 34 ) , A 3 = -id^{v, r 34 ) , A, = id 3 <j>{v, v<®, r 34 ) . (4.3) 

The gauge field configuration A = A^ + A^ has topological charge v = 1 and the 
transition functions are given by 

U 1 = U 3 = 1, U 2 = exp (2™ (g + g) ) > U 4 = exp ^Trig^ . (4.4) 

In this case the eigenvectors of the Dirac Hamiltonian H do not factorize as in the examples 
considered before. In a numerical study we obtained a spectral flow of left handed four- 
spinors from positive to negative energy in accordance with the existence of a left-handed 
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Figure 9: Spectral flow of the left handed eigenmodes of the Dirac Hamiltonian for 
non orthogonally intersecting vortices depending on X4/L4. The vortices cross at 
X4 = L4/2, i.e. at r = r = 0.5. 

zero mode of the Dirac operator in D = 4, see fig. 01 From the figure |H1 one can argue that 
(as in the case of orthogonally intersecting vortices) the largest gradient of the spectral 
flow with respect to time r is where the vortices intersect, r = r , i.e. where the topological 
charge is localized. Further it seems that all eigenmodes cross zero energy level as time 
t evolves from —00 to +00. This is different from the case of orthogonally intersecting 
vortices considered before. There one has different branches in the spectral flow, compare 
eq. ()3. 10p . and only eigenmodes with m = cross zero energy as time r evolves from —00 
to +00. The probability densities of the eigenvectors of the Dirac Hamiltonian again show 
an enhancement at the position of the vortices, see fig. EH 

5 Concluding remarks 

We have studied the quark spectrum in the background of intersecting flat vortex sheets. 
Choosing the four-torus T 4 as space-time manifold we have calculated the spectral flow of 
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Figure 10: Probability density of the eigenmode of the Dirac Hamiltonian whose energy 
eigenvalue crosses zero where the parameters in eq. 1)4.1 Jl have been chosen to be = 

0. 3 + iO.b, = 0.3 at x% = 0.3L 3 , x 4 = 0.29 and 23 = 0.3, 24 = 0.52, respectively. 

the Dirac Hamiltonian as function of time for intersecting vortex fields. The spectral flow 
as function of time r has the largest gradient in the neighborhood of the vortex intersection 
point, i.e. where the topological charge is localized. For the infinitely thick intersecting 
"vortices" considered in section El this gradient is constant (see eq. 1)2.22)) ) reflecting the fact 
that in this case the topological charge is homogeneously distributed over space-time. These 
results are consistent with the findings of ref. JH] that zero modes of the 4-dimensional 
Dirac operator are localized at the intersection points of vortices. Also, the probability 
density of the quark eigenmodes of the Dirac Hamiltonian considered in the present paper 
show a strong correlation to the position of the vortex sheets. Depending on the energy 
level considered, the vortices seem to attract or repel the quarks. For the lowest lying 
states the probability density always shows an enhancement at the position of the vortices, 

1. e. the vortices "act attractively" on the quarks in these modes. 
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A Conventions 

We choose the generators of the gauge group to be anti-hermitian. Therefore the com- 
ponents of the gauge potential are anti-hermitian, e.g. purely imaginary for the gauge 
group U(l). The magnetic flux $ through a closed loop C is defined by 

$= ^i^ (A - i} 

and thus real valued. 
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We consider the Dirac equation in Euclidean space-time. In D = 4 we use the chiral 
representation for the Dirac matrices: 

* = C o)' ! ^ 2 - 3 ' o)' (A - 2) 

75 = 7i727s74 = ( J -1 J ' ( A - 3 ) 

where 1 is the 2x2 unit matrix and Oi are the (hermitian) Pauli matrices. 
The massless Dirac Hamiltonian H in three dimensions reads 

U - - W ,D, . ( ^ ) , (A.4) 

where Dj = di + A4 is the covariant derivative. 

In this paper we use the theta function conventions introduced in |21j : 

9(z, if) = e~^ n2+27rinz , zeC,(G M + . (A.5) 

These functions are analytic in z and obey the periodicity properties 

9(z + 1, iO = iO , 0(* + iO = e^- 2 ™9(z, i£) . (A.6) 
The only zeros of this function are at the points [21] 

z = (m+l/2) + (n + l/2)i£, m,neZ. (A.7) 

B Zero mode of the Dirac operator in D = 4 

We consider the gauge potential 

Ai = -2 7 rix 2 /(L 1 L 2 ) , A 2 = 0, (B.l) 
A 3 = -27rix 4 /(L 3 L A ) , A 4 = (B.2) 

on the four-torus with transition functions 

U 1 = U 3 = 1, U 2 = exp (-27rixi/Li) , C/ 4 = exp (-2mx 3 /L 3 ) . (B.3) 

The topological charge of this configuration is 1, i.e. there is one (left handed) zero mode 
of the corresponding Dirac operator which can be written down explicitly in terms of 
theta functions. To this end one introduces complex variables u := (x\ + ix 2 ) jL\ , v := 
(x 3 + ix&) / L 3 and positive real numbers r 12 := L 2 /Li, r 34 := L 4 /L 3 . The problem of 
finding the zero mode becomes a problem of finding the zero mode of the Dirac operator 
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P(2) = i^iDi + ICT2D2 on the two-torus and this problem can be easily solved (compare 
[To| 120] ). The zero mode of D^ is given by 



\ „ , r r x / ttxI 



,. ^0(^1,^2,^1,-^2) =exp - — — )6(u,iL 2 /L 1 ) . (B.4) 

Vol^l) ^2, -^2j / V -t^l-^2 



A short calculation shows that the zero mode in D = 4 dimensions has the form: 



( \ 




\ <£o(zl,x2,£l, £2)^0(^3,^4, £3, £4) / 



(B.5) 



Obviously, the non-zero spinor component of ipo is the product of the two zero mode 
functions tpo from D = 2. 
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